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We investigate the soliton structure of novel (2+l)-dimensional nonlinear partial differential evolu- 
tion (NLPDE) equations which may govern the behavior of a barothropic relaxing medium beneath 
, high-frequency perturbations. As a result, we may derive some soliton solutions amongst which 

. three typical pattern formations with loop-, cusp- and hump-like shapes. 

o 

^ • In real natural world, there exist very complicated phenomena closely related to nonlinear systems. However, 
in various cases, the real natural phenomena are too intricate to be described only by virtue of (l-l-l)-dimensional 
. NLPDE equations. Thus, a great deal of interest has recently been paid to higher dimensional NLPDE equations, 
' especially, the (2-|-l)-dimensional cases ([ij and references therein). Such higher dimensional systems have been 
, shown to possess soliton solutions which may arise as a balance between nonlinearity and dispersion. One underlying 
query that may arise, consists of the survival of the soliton properties of solutions to such systems with higher order 
O |. nonlinearity. In the wake of such query, we consider a barothropic relaxing medium p — p(p, A) under high-frequency 
' perturbations p' <^ po- The quantities p and p may stand for pressure and mass density, respectively. The constant 
I A is an additional parameter. The pressure po is measured at the unperturbed state. For X = 1, p = pf, inner 
^ ■ interaction processes are frozen and the velocities for fast processes may be defined as v'j = For X — 0, p — Pe, 
I ^ I ■ there is local thermodynamic equilibrium and the velocities for slow processes may be defined as = . 

A few years ago, Danylenko et al. have proposed dynamical equations of state for multicomponent relaxing 
media. In the wake of such result, after a suitable expansion of the specific volume V = as power series of 
, perturbation p' with accuracy o{p'^), the following dynamical equation may be written down 

5 : rdt {dl^^^p' - \^d1pl afd\p'^ -f afd\pl^ + dl^^^pl - \dlp' + a^dlp''' a.dlp'^ =0, (m = 1, • • • , iV), (1) 



, V / 

, where r is the relaxation time; constants and a/ represent positive- valued second order expansion coefficients of 

J Ve and Vf, respectively; constants ae and a/ stand for positive- valued third order expansion coefficients of 14 and 

Q ' V/, respectively; quantity = {x,y,z,- ■■) may stand for position- vector and d^^^,^ = d^^^i + + + ■■ ■ , 
^ {m — 1, ■ ■ ■ , N). It is noted that any repeated index refers to summation with respect to Einstein's notation. In order 
• w^ ', to investigate the eq. H]), the multiscale method [1, |3| may be useful. Thus, defining the quantity e = tlj {lu being 
• the frequency of the processes) chosen to be small (large) parameter, after introducing the following independent 
^ ; variables Tq = tiv, T_4 = twe^'*, X^" = x"'uj, X"_\ = a;"we-'*, (m = I,-- - ,iV) into eq. seven coupled 

equations may be derived. From this coupled system, there may be two leading equations expressed in terms of Tp 
and X™ (m — 1, • • • , N), describing low- frequency perturbations, and the two other expressed in terms of T_4 and 
X!"4 (m = 1, • • • tN) describing high-frequency perturbations. Thus, focusing our interest only to high frequency 
perturbations, we may derive the following evolution equation 

(a,„ + a,„J(a.„ + O^J {p' + afv}p'^ +afv}p'^) - vj^d^p' + I3fd.,,j'^ + -ifp' = 0, {n< m), (2) 

where J™ = p' {ni = 1, • • • , N), and the quantities /3/ and 7/ may be expressed as follows 
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Eq. ([2]) may be obtained in the following way. A dispersion relation for the linearized eq. ^ may be written down 
in the form vj'^u)^ = 

) + jPf J2m=i ^ra " If with (n < m = 1, • • • , TV) . The nonlinear terms may 
be reconstructed in agreement with the initial equation. We may note here that eq. fS]) has dissipative /3/-terms 
and dispersive 7/-terms. If a/ = a/ = /?/ = 0, eq. ^ reduces to a typical modified Klein-Gordon equation in 
(N-t-l)-dimensional space. We particularly focus our interests to = 2 case. Then, it comes 

(9. + dyfp' - vfdlp' + asv}{d^ + dyfp'^ + afv}{d^ + dyfp'^ + I3f{d^ + dy)p' + -ffp' = 0. (4) 

If a/ = and y-terms removed, eq. (jl]) may reduce to the (l+l)-dimcnsional Vakhnenko (V) equation 0, S 0, B 0] ■ 
We need the following accuracy 

(9. + dy)^ - vj^dl « 2(9, + dy) (d^ + dy + vJ^Ot) , (s) 

to further investigate the eq. (jl]). Thus, it may be interesting to consider two important cases: a/ = and a/ 7^ 0. 

1. Case a/ — 0, 

then, eq. (|4]) may be transformed to 



{d, + dy) 

provided 



di-^u^d, + dy) 



u + a{dx + dy)u — u = 0, (6) 



i = - Vft), y = -,J^{y ~ Vft), i^vf^^t, u^ajvjp', a = -^=, (7) 
hold. Performing variable transformation, we introduce new independent variables X, Ti and T2 as follows 

i^Ti-i/" U^dX' + io, y = T2-\l U^dX' + yo, i ^ X, (8) 

where io and j/o stand for arbitrary constants and y, t) — U{Ti, T2, X). Eq. ([6]) may be rewritten as follows 

UxT, + UxT, + &{Un + C/t. ) - (1 + <^ + 0)C/ - 0, (9) 

where (p = — UUxidX' and (f> = — UUT2dX' . Subscripts with respect to X, Ti and T2 may denote 
partial differentiations. Using the ansates ip = Z\Tx ^ 1 ^nd = 2^2T2 ^ Ij cq. ([6]) may bilinearize to 

{DxDt,+DxDt.,+&[Dt,+Bt;)-1)G-F D\F ■ F ^hj^ , (10) 



provided 



C/=^, Zi=Ti-2(ln(F))^, Z2=T2-2(ln(F))^, (11) 



hold. Notations Dx , F)t^ and denote Hirota operators [10, [lH ■ Expanding F and G in a suitable formal 
power series, a one-soliton solution to eq. ([9]) may be given by 

U = 2Ksech{9), Zi = Ti - 2A: [tanh(6l) + 1] , Z2 = r2 - 2Js: [tanh(6l) + 1] . (12) 

where 9 = KX — uJiTi — uJiTi + 9o, 9q being an arbitrary constant. The dispersion relation may be given by 

{K + a){LUi+uj2) + l^0. (13) 

As a result, assuming that uji — Kvi and uj2 — Kv2 such that u = wi + W2 < 0, we may find out three typical 
pattern formations. Indeed, 

(a) for a < J ^^^^ loop-like pattern may be obtained (see FIG. 1 in the case of u = —0.24); 
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(b) for a — etc — y^^: cusp-like pattern may be obtained (see FIG. 2 in the case oi v — —0.24 where 
ttc = 1,4433756729740644112728719512549); 

(c) finally, for a. > \J hump-like pattern may be obtained (see FIG. 3 in the case of v = —0.24). 

These illustrative curves are plotted at initial time i = 0. An essential remark on eq. ([6]) should be noted here. 
In the wake of the results got from ref. [12], an extended complex- valued form of eq. ^ may be found as follows 



{di + dy) 



q + a{dj. + dy)q'- g = 0, 



(14) 



where q ^ + i<z™ may stand for a complex-valued observable, symbol (★) may refer to complex conjugation 
and = —1. Thus, one may derive the following set 



(15a) 



QxT, + QTt. + 5 (OtT + Qri) - (1 + ZiT, + Z^T.)Q"'^ = 0, 



(15b) 



where q(i, y, t) = Q{Ti, T2, X) such that independent variables x, y and t now stand for 

S: = Ti-\j QQ*dX' + io, y = T2 - i / QQ*dX' + yo, i=X. 



(15c) 



(16) 



It comes that eq. maybe closely related to the coupled dispersionless systems [H, [T3] recently investigated 
by Kakuhata and Konno [iHllBl- Besides, eq. ([TS]) may be also observed as a coupled (2-|-l)-dimensional version 
of the complex-form of the Schdfer- Wayne short pulse (SWSP) equation 17] that has been subject to many 
recent investigations [13, E, 113, HII ■ A one-soliton solution to eq. (fT5|) may be written as follows 



Q = Aexp(6'™)sec/i(6l''), = Ti - 2K' [tanh(0'') + 1] , ^2 = 72 - [tanh(6l'" 



(17) 



where 6 — 6"^ + lO™'^ = KX — TiUJi — T2UJ2 + ^Oi ^'o being an arbitrary complex-valued constant, K = K^' + iK^ 



and uji 



[j — 1,2). The amplitude A may be given hy A = 2K^. The corresponding dispersion 



relation may be given by eq. (|13p . As a result, one may easily find that solutions given by eq. (|17p in terms of 
and Q™, possess a nonzero angular momentum i.e. the previous patterns depicted above may rotate with 
an angular frequency SI — ojf"' -I- 1^^™- Concretely, if we assume that — viK^ , K^™ = C'^i™, — V2K^, 
j^im _ i^Qj]^ V — vi + V2 < 0, and Q. may be expressed as follows 



vcV (c--i)^' 



(18) 



provided < C < Thus, setting q, 



— — ^ and (5s — t^-tt^, it comes that, 



/2\v\ 



(a) for a > (5s, no rotating pattern formation is expected; 

(b) for oic < a < (5s, rotating loop-like pattern formation is derived (see FIG.l for v = —0.24 and ( = 13/3); 

(c) for a = Uc, rotating cusp-like pattern formation is derived (see FIG. 2 for v = —0.24 and ( — 13/3 where 
ttc = 0,057735026918962576450914878050196); 

(d) finally, for a < etc, rotating hump-like pattern formation is derived (see FIG. 3 for v = —0.24 and C = 13/3). 

2. Case a/ ^ 0, 

then, eq. (HI) may lead to 



(di + dy) 



di + u{dx + dy) + -u {di + dy) 



u + a{dx + dy)u + u = 0, 



(19) 
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provided 

^ = y = \^miy-^f*)^ i^\f7^^f^}t^ ^ = —p'^ fi = ^W^, (20) 

Y M^/ V f f V 6a/ ^ af afVf \j 27/ 

hold. Wc introduce new independent variables X', T[ and as follows 

~x = T[ + \( {U^ + 2U)dS + xi y = n + \l {U^ + 2U)dS + yi i^X', (21) 

where x'q and j/g stand for arbitrary constants and u(x,y,t) = U(T{,T2, X'). Then eq. (fT9|) becomes 

Ux'Ti + Ux'Ti + a{UTi + Uri) + {1 + ^' + ^')U = 0, (22) 

where if' = J^^{U+l)UT,dS and 0' = /_^^(C/+l)C/T,d5'. Using the ansates (p' = -Zit[ + ^ and 0' = -Zar^ + l, 
eq. (Uni) may be bilinearized as 

{Dx'DT^+Dx'DT^+a{DT^+DT^) + l)G-F^Q, D^F ■ F = ^ {G^ + 2GF) , (23) 

provided 

U^^, Zi=r{-2(ln(F))^,, Z2 = r^-2(ln(F))^,, (24) 

hold. Thus expanding suitably the functions G and F in power series according to Hirota method [13, El, 
one-soliton solutions of kink-like pattern may be expected and discussed by means of the dispersion relation 
derived from eq. ((23|l . 

In summary, we have given a general model of a (2+l)-dimensional NLPDE equation (see eq. (fT4|) ) with the soliton 
structure. This novel (2-|-l)-dimensional NLPDE equation (eq. (fT4|) ) may be observed as a coupled (2+l)-dimensional 
version of the (l-l-l)-dimensional complex-SWSP equation This novel equation with the other one (see eq. 
may deserve further scientific interests both from the viewpoint of the investigation of the propagation of high- 
frequency perturbations and from the viewpoint of the existence of stable wave formations. We may actually believe 
that these equations may be valuable for studies on soliton theory, geodynamics, plasma physics, hydrodynamics, 
condensed matter, string theory, nonlinear optics, just to name a few. 
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